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Abstract. We present early WHT ISIS optical spectroscopy of the afterglow of gamma-ray burst GRB 050730. The spectrum 
shows a DLA system with the highest measured hydrogen column to date: N(R I) = 22. 1 ± 0. 1 at the third-highest GRB redshift 
z = 3.968. Our analysis of the SwiftXKF X-ray observations of the early afterglow show X-ray flares accompanied by decreasing 
X-ray absorption. From both the optical and the X-ray spectra we constrain the dust and gas properties of the host galaxy. We 
find the host to be a low metallicity galaxy, with low dust content. Much of the X-ray absorbing gas is situated close to the 
GRB, whilst the H I absorption causing the DLA is most likely located further out. 
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1. Introduction 

Gamma-ray bursts (GRBs) have proven to be excellent probes 
of the distant Universe. High luminosity GRB afterglows allow 
absorption line studies of the ISM at high redshift to at least 
Z = 4.5 (see Andersen et al. 2000). The launch and success- 
ful operation of the Swift satellite means more GRBs are be- 
ing localised and afterglows studied. Subsequently, the number 
of high redshift bursts suitable for host galaxy spectral stud- 
ies has dramatically increased. Deep observations of afterglow 
positions have detected host galaxies in almost all cases (e.g. 
Conselice et al. 2005). Most hosts are compact, actively star- 
forming galaxies and, where the relevant data are available, are 
found to have low metallicity and low intrinsic extinction (e.g. 
Berger et al. 2003; Tanvir et al. 2004; Christensen et al. 2004). 
However, in a few cases, radio/submm observations of hosts 
give a star-formation rate (SFR) which is of order a few to 
~100 times larger than rates derived from optical estimators 
such as the line luminosities of Ha and [O II] or the 2800 A 
restframe UV continuum flux (e.g. Berger et al. 2003). This 
may be caused by strong dust obscuration, but neither spectra 
nor colours of hosts show strong internal extinction. Afterglow 
spectroscopy provides a unique window on the near environ- 
ment of GRBs (e.g. GRB 021004, Schaefer et al. 2003; Fiore 



et al. 2005; Starling et al. 2005), allowing us to probe the ab- 
sorbing dust and gas properties in more detail. In this Letter 
we present optical and X-ray spectra of GRB 050730, discov- 
ered by Swift on July 30th 2005, 19:58:23 UT (Holland et al. 
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2005) and lying at a redshift of z = 3.97 (Chen et al. 2005] Rol 
et al. 2005), in which we study the circumburst gas and dust 
properties. 

2. The optical afterglow spectra 

2.1. Observations 

During the afterglow phase of GRB 050730, we acquired 
spectra using the Intermediate-dispersion Spectroscopic and 
Imaging System (ISIS) on the William Herschel Telescope. 
The R316R and R300B grisms were used on the red and blue 
arms respectively. Two observations were done sequentially, at 
the parallactic angle, with exposure times of 1260 and 1800 
seconds. The first observation started at 22:57 UT at airmass 
~2.73 (midpoint 0.132 days after burst), the second at 23:19 
UT (midpoint 0.145 d) and airmass ~3.4. The seeing quality 
at the high airmasses required the slit width to be widened to 
2.5 arcsec. Conditions during the observations were not photo- 
metric. These factors mean that our absolute flux calibration is 
not reliable, but the relative calibration should not be affected. 
Both spectra have been reduced using the data reduction pack- 
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Fig. 1. The WHT ISIS combined, normalised spectrum of the afterglow (midpt 0. 14 days), and lcr error spectrum (lower curve). 
Overlaid is the best-fitting DLA profile (solid line) and its errors (dashed lines). All significant lines (3cr) are indicated above 
for z = 3.969 (solid), z = 3.565 (dashed), z = 1 -773 (dot-dashed) and unidentified (dotted) systems; see on-line table for details. 
Lines used in further analysis are labelled. 



age IRAF following standard procedures. A Galactic extinction 
correction of E(B - V) = 0.049 (Schlegel et al. [T9^St was ap- 
plied. The wavelength resolutions of blue- and red-arm spectra 
respectively are 8.7 and 8.1 A. The signal to noise per pixel, 
measured at 6800 A, is 27 in the first and 17 in the second 
spectrum. 

2.2. Results 

The spectrum, shown in Fig. ^ is rich i n li ne features at 
z = 3.97,3.56 and 1.77. A strong Damped Lyman-Alpha ab- 
sorption system (DLA) is present; here we focus on this and a 
selection of metal lines presumed to originate in the GRB host 
galaxy. We fitted a power law continuum corrected for Galactic 
extinction to the ~6500-7500 A region of each spectrum, ex- 
cluding the absorption lines, and find an epoch averaged slope 
of B — -1.34 + 0.21 (2cr formal fit error). We tested for any 
departure from a pure power law due to host-galaxy extinction: 
fitting MW, LMC and SMC extinction curves (Pei 1992) all re- 
sult in epoch averaged Ay = 0.01. The optical/IR spectral slope 
from published BVRIJ photometry extrapolated to a common 
epoch using a temporal decay slope of 0.89 (Haislip et al. 2005 
Holman et al. 120051 Cobb et al. l2005l Blustin et al. l2005> gives 
B 1 , consistent with the spectral analysis but not very con- 
straining. 

Despite the moderate dispersion of the ISIS 300 grisms, 
the damping wings of the host galaxy DLA are clearly visi- 
ble. In fact, the determination of N(H I) in DLAs based on long 
slit spectra is considerably simpler than for echelles. Since the 
damped profile may extend over many spectral orders in a typ- 
ical echelle, accurate combination and flux calibration can be 
troublesome. Using the Starlink software DIPSO, we determine 
log AT(HI) = 22.1 +0.1 (see Fig. [TJ. Taking LyB into account 
did not lead to a more accurate determination of N(H I), and the 
error on our fit is dominated by uncertainties in the determina- 



tion of the power law continuum. The A^(H I) value is consis- 
tent with, although slightly lower than, that reported by Chen 
et al. (2005). This high value (the first DLA to break the 10 22 
atoms cirT 2 barrier) continues the trend amongst GRB DLAs 
towards very high neutral hydrogen columns (e.g. Jensen et al. 
2001; Hjorth et al. 2003; Vreeswijk et al. 2004). 

Although our spectra do not enable as detailed a study of 
the metal lines as is possible via echelle observations (e.g. 
Chen et al. 2005), we briefly comment on a selection of these. 
Detection limits are quoted at the 3 cr level. Although we detect 
both S II A, A 1253, 1259, both lines are likely to be at least par- 
tially saturated. In addition, the weaker S II A 1253 line which 
potentially offers a better limit on A^(S II) is blended with an- 
other (unidentified) feature (H.-W. Chen, 2005, private commu- 
nication). We determine an upper limit of [S/H] < -2.0 based 
on the absence of the weaker S II A 1250 A line 1 , in good agree- 
ment with Chen et al. (2005). Similarly, from the Fell /U608 
line which is partially saturated and the undetected Fe II A161 1, 
we determine -2.9 < [Fe/H] < -1.9. 

A search for variability in line features between our two spec- 
tra, separated by ~25 mins, revealed no significant changes (see 
on-line table for details), neither did the N(H I) column vary. 



3. The absorbed X-ray afterglow 

3.1. Observations 

We have analysed the early Swift XRT data, to look for evi- 
dence of intrinsic absorption in the X-ray spectrum. The XRT 
data consist of Windowed Timing (WT) mode data for the first 
orbit (133 to 793 seconds after the trigger) and start of the sec- 
ond orbit, and Photon Counting (PC) mode data for later orbits. 
The data were reduced using the standard pipeline for XRT 
data within the HEADAS 6.0 package (Swift software version 



S abundance of log (S/H) + 12 = 7.20 (Grevesse & Sauval 1998) 
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Fig. 2. Evolution of the 0.3-10 keV count rate (top), 0.2-10 
keV power law slope, T, (middle, note that r = 1 -J3) and addi- 
tional equivalent hydrogen column, Nmnt, at the redshift of the 
host galaxy (bottom) during the first ~800 s of the Swift XRT 
observations (90% errors). The individual spectra are shown 
with filled circles; the two combined spectra pre- and post-500 
s, for which better constraints are obtained, are shown with di- 
amonds. 



2.0). WT mode data was extracted using a rectangular region 
centred on the source, and a similar area in a source-free re- 
gion of the same image to determine the background level. PC 
mode data was extracted using a circular aperture, except for 
orbits 2 to 4 which show evidence of pile-up (count rates > 0.8 
counts s _1 ) and were extracted using an annular region cen- 
tred on the source and filtered on grade only. The light curve 
was obtained between channels 30 and 1000 (spanning ~0.3- 
10 keV). Spectral analysis was done using XSpec 11.3, with the 
standard Ancillary Response Function (ARF) files, which esti- 
mate the effective telescope area, for PC mode data, and with 
ARFs based on ray-tracing ('physical' ARFs) for WT mode 
data which should provide a better calibration at low energies. 

3.2. Results 

The first orbit shows several flares in the light curve, first re- 
ported for this afterglow by Grupe et al. (2005). We have 
performed a detailed analysis of the spectral evolution of the 
early-time data. The fitted model consists of a power law plus 
Galactic absorption (fixed at 3.05 x 10 20 cm~ 2 , Dickey & 
Lockman ll990> and a variable Galactic-like absorption com- 
ponent with Solar metallicity and z = 3.97. Errors are quoted 
at the 90% confidence level for 1 interesting parameter. We 
find evidence for a change in power law photon index, from 
r = 1.52 ± 0.04 at the start of the first orbit to 1 .79 + 0.06 at the 
end of the orbit (note that r = 1—0). We also find evidence for 
an excess absorption column, which at the redshift of the burst 
amounts to an intrinsic column of Nhm = (1.4 ± 0.3) x 10 22 



cm . However, around 500 seconds post trigger, the absorp- 
tion column abruptly changes, becoming lower by about a fac- 
tor of 4: A^H.int = (3.4 + 2.7) x 10 21 cm" 2 . In the late-time 
PC mode spectrum the intrinsic column cannot be constrained, 
setting an upper limit of only Nu,m < 1.0 X 10 22 crrT 2 , and 
the power law photon index remains stable at T ~ 1.77. We 
have checked for a possible correlation between the intrinsic 
Nn.i n t and T in the fit. Contour plots for the intervals 133-503 s 
and 503-793 s post trigger show no evidence for any correla- 
tion, confirming the reality of both the drop in Nhm and in- 
crease in T (Fig. E}. Interestingly, this happens directly after 
the peak of the second visible flare, where the light curve in- 
tensity has increased by a factor of 3. Given the host galaxy 
metallicity we measure in the optical spectrum, we adjust the 
X-ray absorption model accordingly. Using Z = Z o /100 for 
all the elements heavier than He included in the zvphabs X-ray 
absorption model, the required intrinsic equivalent hydrogen 
column increases by a factor of ~10 in both cases to A^int = 
9.5! 2 ] x 10 22 cnr 2 (first 400 s) and N HM = 2.6+Jj? x 10 22 cnT 2 
(> 500 s post trigger) with approximately the same goodness 
of fit. 

A preliminary analysis of published optical photometry to- 
gether with the PC mode XRT spectrum has shown the X-ray 



and optical slopes at 0.19 days (J3 X ~ -0.7 to -0.8,/?, 



opt 



-1.0 



to -1.5) to likely be incompatible with a position of the cool- 
ing break between the optical and xrays, and might suggest the 
presence of an inverse Compton component; we await the avail- 
ability of further optical/IR photometry for a full analysis. 

4. Discussion and conclusions 

4. 1 . Host galaxy properties 

There is a well known relationship between galaxy luminos- 
ity and metallicity (e.g. Garnett 2002; Lamareille et al. 2004) 
which spans 6 orders or more of magnitude in Mb- Tremonti 
et al. (2004) have recently demonstrated that this relation is 
driven by an underlying relation between mass and metallic- 
ity. The cause of the relationship, they argue, is due to the in- 
creased gravitational potential of massive galaxies which en- 
hances metal retention. In the absence of a detected host for 
GRB 050730 at the time of writing, it is in principle possi- 
ble to use the luminosity-metallicity (LZ) and mass-metallicity 
(MZ) relations to predict the Mb and stellar mass of the 
host. Both of these relations are best determined locally (e.g. 
Lamareille et al. 2004), although sizeable datasets have now 
investigated the LZ relation up to z ~ 1 (e.g. Kobulnicky et 
al. 2003; Kobulnicky & Kewley 2004). There is clear evidence 
for evolution in the LZ relation, in the sense that galaxies are 
more metal-poor for their luminosity at higher z (although see 
caveats in Kewley & Ellison in prep.). This trend appears to 
continue both for the LZ and MZ relations up to z ~ 3 (e.g. 
Shapley et al. 2004; M0ller et al. 2004; Erb et al. in prep.), 
although only the highest mass/luminosity galaxies are bright 
enough to be included in spectroscopic samples. The lowest 
metallicity bin in the fitted MZ relation of Erb et al. (in prep.) 
is Z ~ Z /3 corresponding to a stellar mass log(M*/M ) ~ 9.5. 
The metallicity measured from absorption lines in the optical 
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afterglow considered here is Z ~ Z G /100, which indicates that 
the host is not a massive, luminous Lyman break galaxy (LBG), 
although Jakobsson et al. (2005) argue that GRB hosts follow 
the same UV luminosity function as the faint LBGs. We do note 
that the MZ relation is based on emission lines. However, HST 
imaging has shown that GRBs occur in regions of strongest 
star formation (e.g. Fruchter et al. 2005), justifying our as- 
sumption that the absorption lines are formed in the same re- 
gions as the higher wavelength emission lines. Combining the 
measured N(H I) with the metallicity and assuming an SMC 
gas-to-reddening ratio (Bouchet et al. 1985), we can estimate 
the extinction associated with the GRB host galaxy and com- 
pare this to the values obtained from the optical continuum fits. 
Negligible E(B - V) is determined by both methods, consistent 
with the similarly small amounts of dust seen towards inter- 
vening DLAs (e.g. Murphy & Liske 2004; Ellison, Hall & Lira 
2005). 

4.2. The neutral hydrogen column 

GRB 050730 has the strongest DLA seen in a GRB afterglow 
spectrum, with a hydrogen column density of log N(H I) = 22.1 
+ 0.1. The X-ray absorption at late times scaled to Z o /100 
yields a comparable log A^im = 22.4^^ (assuming the Nn,mt 
measured at -500-800 s post burst in the WT mode XRT spec- 
trum can be extrapolated to a few hours post burst). The Nhm 
we measure in the early-time X-ray spectra covering ~133- 
500 s post trigger is about ten times higher than that measured 
at t > 500 s. The change in X-ray absorbing column could be 
caused by ionisation by the gamma-ray jet, or by the X-ray 
flares which are suggested to be caused by prolonged central 
engine activity (Burrows et al. 2005; King et al. 2005). 

It should be noted that what is measured in the X-ray mod- 
els is an equivalent hydrogen column, since primarily metal 
edges contribute to the X-ray absorption at the redshift of 
GRB 050730, and that this is highly dependent upon the metal- 
licity assumed (see e.g. Wilms et al. 2000). There will be a 
contribution to the X-ray absorption from intervening systems, 
which cannot be disentangled from absorption in the host, par- 
ticularly given that we do not know the metallicity of the clos- 
est intervening system observed in this spectrum (z = 1.77). 
In principle, a lower column very close to the observer could 
have a similar effect on the spectrum as a large column at high 
redshift. 

The observed X-ray column variability does, however, 
lead us to conclude that most of the X-ray absorbing gas in 
GRB 050730 is located close to the GRB. The optical H I col- 
umn remained stable over the ~25 mins between our ISIS spec- 
tra, taken at 0.132 days since burst, well after the observed 
X-ray flaring (although the occurrence of X-ray flares at later 
times cannot be ruled out owing to low count rates). The H I 
creating the DLA is likely to be located much further away 
from the GRB, unaffected by the GRB radiation. We would 
expect to observe destruction by the GRB of dust co-located 
with the X-ray absorbing gas. Our spectra imply a very low ex- 
tinction in the host at ~3 hrs post burst. Future prompt optical 



Table 1. Lines detected above 3<x in the first and second epoch 
WHT spectra. 3 . Submitted as on-line table only. 
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0. 
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09 
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0. 


12 
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3 


.5641 
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a Blueward of Lyo- the low resolution and Lye forest hampers se- 
cure identification of metal lines, which we therefore do not list. 

b Due to the uncertain continuum level, we do not attempt to mea- 
sure the widths of the Lyo- and Ly/J lines. 

c This equivalent width measurement is seriously affected or made 
impossible by the atmospheric absorption band from 7584- 
7675A. 



spectra, in conjunction with X-ray observations, are required to 
investigate this further. 
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